Cuts and penalties: a comment on "The clustering of ultra-high energy cosmic rays 

and their sources" 



O 

o 

(N 



(N 
> 

m 
m 

o 
m 
o 

Oh 
I 

o ■ 



P.G. Tinyakov'^-'^ and I.I. Tkachev*'''= 
""Institute of Theoretical Physics, University of Lausanne, 

CH-1015 Lausanne, Switzerland 
^CERN Theory Division, CH-1211 Geneva 23, Switzerland 
'^Institute for Nuclear Research, Moscow 117312, Russia 

In a series of papers we have found statistically significant correlations between arrival directions 
of ultra-high energy cosmic rays and BL Lacertae objects. Recently, our calculations were partly 
repeated by Evans, Ferrer and Sarkar Q with different conclusions. We demonstrate that the 
criticism of Ref. Q] is incorrect. We also present the details of our method. 



I. INTRODUCTION 

Identification of sources of the ultra-high energy cosmic 
rays (UHECR) is a key to understanding their nature. 
The observed smalLscale clustering of UHECR 0,111113 
suggests that already existing data 0,0, contain 
information sufficient to make first steps in this direc- 
tion. In a series of papers jSli|l3,|ll|ll[ll|ll we 
have shown that there exist significant correlations be- 
tween the arrival directions of UHECR and BL Lacertae 
objects (BL Lacs), and therefore BL Lacs are likely to 
be sources of UHECR. Although present small dataset 
does not allow to determine with certainty the details of 
UHECR production and propagation, interesting hints 
may be obtained [H El 111 111 . 

The first evidence of connection between UHECR and 
BL Lacs was found in Ref. @| where we have shown 
that the combined set of AGASA events with E > 
4.8 X 10^^ eV and Yakutsk events with E>2Ax 10^^ eV 
(identified previously as the sets with largest clustering) 
correlates strongly with most powerful BL Lacs. The 
statistical significance of this correlation was found to be 
6 X 10~^ with penalty factor included (for explanation of 
penalty factor see below). 

In the recent paper by Evans, Ferrer and Sarkar (EES) 
the validity of this result was called into question. The 
rational part of the criticism of EFS boils down to the 
following two issues: i) choice of cuts on BL Lacs and/or 
calculation of the penalty factor ii) choice of the cosmic 
ray set. Below we address these issues and demonstrate 
that criticism of EFS is unjustified. We also clarify some 
frequently arising questions related to our analysis. 



II. STATISTICAL SIGNIFICANCE OF 
CORRELATIONS: CUTS AND PENALTY 
FACTORS 

A. What is penalty factor 

In statistical analysis, one nearly always has to make 
cuts in order to improve signal-to-noise ratio. The ques- 
tion is how to take them into account correctly. A com- 



mon wisdom is that cuts should be fixed a priori, i.e., 
based on theoretical considerations. In that case the cuts 
simply limit the number of data points, but do not alter 
the calculation of probabilities. 

Not always this is possible. For instance, in the case 
of UHECRs, their acceleration mechanism is not known. 
How to impose cuts in a catalogue of astrophysical ob- 
jects in order to select actual UHECR emitters? In which 
band — radio, optical. X-ray, 7, TeV? To which flux 
limit? In Ref. [g we proposed an approach to this prob- 
lem which consists in adjusting cut(s) so as to maximize 
the signal, and then compensating this cut adjustment 
by a penalty factor. 

A penalty factor takes into account the "number of in- 
dependent attempts" made when searching for the best 
signal. For instance, if two independent catalogues were 
tried, the penalty factor is 2 (in the limit of small proba- 
bilities). If N catalogues tried are not independent, as in 
the case of cuts in a single catalogue, the penalty factor 
is smaller than N and should be calculated by means of 
the Monte-Carlo simulation. 

In correlation analysis of cosmic rays, the quantity of 
interest is the probability p'^f^^ that the observed excess 
of cosmic rays around source positions is the result of 
a chance coincidence. This probability depends on cuts 
made in the source catalogue; the latter are adjusted so 
that the probability is minimum (correlations are max- 
imum). The cut adjustment should be compensated by 
the penalty factor which is calculated as follows. A ran- 
dom set of cosmic rays is generated and treated exactly 
as the real data: the same cuts are tried and the resulting 
minimum probability p^;'^ is determined. This procedure 
is then repeated for a large number of random "cosmic 
ray" sets, and p^^ is determined each time. The num- 
ber of occurrences of a value p^^ < P is then counted 
as a function of p. Divided by the total number of sets, 
this gives the probability P{p) that the adjustments of 
cuts produces Pmin < P for a random set of cosmic rays. 
P{p'^i^) is the correct measure of the significance of ob- 
served correlations. 

When no adjustment of cuts is made one obviously has 
Pip) — Pi the quantitative measure of correlation is 
the probability p itself. When cuts are adjusted, small 
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probabilities appear more often by construction, and this 
relation is modified. The modification can be written in 
terms of the penalty factor F{p), 

P{p)^F{p)-p. (1) 

In the presence of cut adjustment, the significance of cor- 
relations is determined by the product -flp^jf^f ) • Pnfhf ■ 

To summarize, all parameters characterizing UHECR 
sources and UHECRs themselves can be subject to cuts. 
Some of these cuts can be decided a priori; they reflect 
(and depend on) physical assumptions made. These cuts 
do not require penalty. One can call these cuts fixed. 
Alternatively, in the absence of physical arguments, the 
cut should be chosen (adjusted) so as to maximize the 
signal. These cuts can be called adjustable. These are the 
adjustable cuts which imply non-trivial penalty factors; 
the latter can be calculated in a way explained above. 



B. UHECR set 

In our analysis Q, the cosmic ray set was fixed on 
physical grounds using results of previous publications 
1^, and was not adjusted in search for correlations. The 
motivation was as follows. It was observed earlier |^|^|^ 
that highest-energy cosmic rays exhibit remarkable prop- 
erty: they auto-correlate on the angular scale consistent 
with the angular resolution of detectors. The mere ex- 
istence of these correlations suggests that a relatively 
small number of point sources of cosmic rays is con- 
tributing a considerable fraction of the UHECR flux , 
and that the trace of the original directions to sources is 
not lost completely. The identification of these sources 
via cross-correlation analysis may therefore be possible. 
Within the hypothesis that UHECR clusters are due to 
point sources, the strongest correlation signal is expected 
with the UHECR set which has most significant auto- 
correlations. The latter requirement selects |^] AGASA 
events with E > 4.8 x 10^^ eV and Yakutsk events with 
E > 2.4 X 10^^ eV. The observed angular size of clusters 
also fixes the angular scale at which correlations with 
sources are expected. 

The fact that AGASA and Yakutsk sets have differ- 
ent energy cuts in this approach is, of course, disturbing. 
It may, however, be explained by different energy cali- 
bration of the two experiments, or smaller number of the 
Yakutsk events. In either case, within our assumptions it 
is inconsistent to change energy cuts or to drop Yakutsk 
events once autocorrelations are found in this set. 

The authors of Ref. 1] have rejected the Yakutsk 
events because Yakutsk array has angular resolution 
not as good as AGASA. By itself, worse angular res- 
olution does not imply that correlations with sources 
must be absent in Yakutsk set: even though the an- 
gular resolution is worse, the density of UHECR events 
around actual sources is larger as compared to a ran- 
dom set, and one has an excess in counts even at 



small angles. The actual presence of autocorrela- 
tions in Yakutsk data supports this statement. There- 
fore, it is not correct that "the correlations found at 
smaller (than Yakutsk angular resolution) angles cannot 
be meaningful" 



C. Correlations with BL Lacs 

The authors of Ref. Q have essentially reproduced our 
Monte-Carlo simulations: solid curve of Fig. 3 of Ref. 
agrees with our calculations if Yakutsk events are dis- 
carded. They concluded, however, that we have under- 
estimated the penalty factor. This conclusion was based 
on i) a comparison of correlations in the case when one 
particular set of cuts was chosen to the case when no cuts 
on BL Lacs were made at all ii) smaller correlations of 
BL Lacs with a different set of UHECR. Neither such a 
comparison, nor the dependence on the UHECR set have 
anything to do with the penalty factor. The UHECR set 
was fixed as described above and was not adjusted during 
the calculation, thus requiring no penalty. The penalty 
associated with the adjustment of BL Lac set should be 
calculated as described above. This calculation was per- 
formed in our paper j^] as follows. 

The most complete catalog of QSO contains 306 
confirmed BL Lacs for which the apparent magnitude, 
redshift (where known) and 6 cm radio flux are listed 
(for several objects the 11 cm radio flux is also given). 
In Ref. S] we considered BL Lacs with z > 0.1 as sug- 
gested by statistics of clustering. (We also examined a 
complimetary small set without further cut adjustment. 
This "independent attempt" adds one to the penalty fac- 
tor.) We adjusted cuts on magnitude and 6 cm radio flux 
in the following way. First, the grid of 11 x 11 cuts was 
fixed. Namely, the cut on magnitude was varied from 
m < 16 to m < 20 with the step 0.4; simultaneously 
the cut on 6 cm radio flux was changed from 0.21 Jy to 
0.01 Jy with the step 0.02 Jy. Note that further "refin- 
ing" of this set of cuts would not change the results as 
the sets of BL Lacs obtained would overlap almost com- 
jletely. The lowest probability (best correlation) found 
^in the real data in this way equals to pjjfj'f = 4 x 10^^ 

To calculate the penalty factor associated with cut ad- 
justment, a random configuration of 65 UHECR events 
was generated. Then its correlation with all 11 x 11 sub- 
sets of BL Lacs was examined and the minimum prob- 
ability Prnin ovcr 11x11 cases was found and recorded. 
The whole procedure was repeated for the next random 
configuration. A total of A'tot — 10^ random configu- 
rations were treated in this way. The number of occur- 
rences N{p) of Pmin < P was Calculated as a function of p. 
Divided by iVtot it determines the corrected probability 
P{p) and the penalty factor F{p) as defined in Eq. 
The resulting penalty factor F{p) is shown in Fig.^ The 
error bars correspond to statistical errors in determina- 
tion of N{p). 
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FIG. 1: MC calculation of the penalty factor Pijp)- Solid line 
represents the fit by a power law in the log-log scale. 

Taking the penalty factor at the minimum probabil- 
ity obtained for the real data, F(4 x 10^®) ~ 14 and 
adding 1 as penalty for BL Lacs with z < 0.1 as ex- 
plained above, one finds the significance of correlations 
including penalty, 

P(4 X 10"^) ~ 6 X 10"^ 

The significance is ~ 4tT in terms of Gaussian distribu- 
tion. It should be pointed out that, in order to get the 
correct significance, the penalty factor should be multi- 
plied by the minimum probability obtained for the real 
data with the same set of cuts. 

We stress again that the absence of correlations of 
cosmic rays with the entire sample of known BL Lacs 
does not mean anything. Physically, cosmic rays corre- 
late with the brightest BL Lacs. Including dimmer BL 
Lacs in the analysis dilutes the correlations, and this is 
the only conclusion one draws from the analysis made in 
Ref. 0. 

D. "Correlations" with GRB 

The authors of Ref. [J have chosen GRBs as a con- 
trol set where correlations are not expected; it has to 
be compared with the case of BL Lacs. From the fact 
that correlations with the whole catalogue are absent in 
both cases while there are many coincidences, the authors 
of Ref. concluded that GRBs correlate with UHECR 
"just as well as do BL Lacs!" , implying that correlations 
with BL Lacs are due to cut adjustment. This conclu- 
sion can be tested by Montc-Carlo simulation: one has 
to check whether cut adjustment leads to apparent corre- 
lations with GRBs. It does not, as we now demonstrate. 

The Fourth BATSE Gamma-Ray Burst Catalog ITtI 
contains 1637 objects from which we cut out, in sequen- 
tial order (according to observation date and time), 5 
non-overlapping sets of 306 objects each (recall that 306 



3 




J L 



306 612 918 1224 1530 



FIG. 2: Correlations of UHECR with GRB: probability of the 
observed excess of UHECR events near GRBs for 5 indepen- 
dent GRB subsets as a function of the cut applied within each 
subset. 



is the number of confirmed BL Lacs in the catalogue p^). 
In each set we have looked for maximum correlations by 
optimizing cuts. We present two different tests: one- 
dimensional cut which allows transparent graphical rep- 
resentation, and two-dimensional cut which mimics more 
precisely the case of BL Lacs. We took the same maxi- 
mally auto-correlated UHECR data set as in Refs. 0,0. 

In the first test we replace the 11 x 11 cuts on mag- 
nitude and radio flux by 121 cut on the total number of 
objects (we take them in sequential order, so that first 
few cuts correspond to including 3, 6, 9 ... first objects, 
and the last cut corresponds to including all 306 objects) 
pl|. Like in the case of BL Lacs, we then calculate the 
probability of the observed excess of UHECRs around 
GRBs as a function of the cut. The results are pre- 
sented in Fig. 121 The smallest probability found in all 
5 sets is > 4% (which by itself is not statistically signifi- 
cant). In view of the similar number of cuts in two cases. 
Fig. 1 provides an estimate for the penalty factor, 
F(0.04) « 5, in a single set of 306 GRBs. Multiplying 
the best probability, the penalty factor and the number 
of independent GRB sets one gets a number of order 1. 

To mimic the case of BL Lacs we perform a similar test 
with two independent cuts. The first cut is done as before 
on the sequential number of the event. The second cut is 
imposed on the time (in seconds of the day, UT) of the 
event. Note that these two cuts are virtually uncorrelated 
since BATSE registered roughly one event per day. We 
have searched for the best correlation signal on the grid 
of 11 X 11 equidistant cuts from 1 to 306 and from 1 to 
86400, respectively. In the above five catalogs consisting 
of 306 GRBs each, the following lowest probabilities were 
found: 9.5%, 1.4%, 20%, 1.8%, 24%. Multiplying the 
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lowest of these probabilities 0.014 by the penalty factor 
F(O.Ol) « 8 and by the number of independent GRB 
sets, we again obtain a number of order 1. 

These examples are a good illustration of compensa- 
tion between the effects due to cut adjustment and the 
penalty factor. When compensated by the penalty fac- 
tor, cut selection does not introduce apparent correla- 
tions when in reality they are absent. 

To avoid confusion, we note again that we used GRBs 
here as a control set, and imposed physically unmotivated 
cuts on purpose. 



III. SUBTLETIES OF CORRELATION 
ANALYSIS 

In this section we discuss some subtleties of correla- 
tion analysis which were not mentioned in Ref. We 
hope that this will answer a number frequently arising 
questions regarding our work and correlation analysis in 
general. 



A. Why completeness of the BL Lac catalogue is 
not necessary for establishing correlations with 
UHECR 

A catalogue of astrophysical objects is called complete 
when it is believed to contain more than, say, 90% of 
all existing objects of a given type, in a given region of 
the sky, and down to some fixed level of luminosity in 
a given waveband(s). The catalogue 0|, used in our 
analysis, is not complete simply because it contains all 
objects known up to date without further selection [l^ . 

It is a common wisdom that completeness of the cat- 
alogue is crucial for statistical analysis. The reason is 
obvious: in an incomplete catalogue, the distribution of 
objects in brightness, position, etc., reflects not only their 
actual abundance, but also observational bias. When 
one studies, for instance, the evolution of abundance 
of some objects with time, one has to make sure that 
logN/logS dependence reflects actual change in spatial 
density rather than the fact that remote (and therefore 
dim) objects are more difficult to observe. Indeed, as 
it is, the catalogue is not suitable for any analysis 
where statistical properties of BL Lacs are of interest. 

In the correlation analysis of UHECR the question is 
different: given a set of candidate sources (BL Lacs in 
the case at hand), is the distribution of cosmic rays ran- 
dom, or does it peak around BL Lac positions? Clearly, 
these are the statistical properties of cosmic rays which 
are of interest in that case, not the statistical properties 
of BL Lacs. Technically speaking, when calculating the 
correlation function by the algorithm of Ref. , the cos- 
mic ray directions are simulated, while source positions 
are held fixed. For this reason the method is applicable 
without change to studying, for instance, a correlation of 
UHECR with one particular direction (say, the Galactic 



center), i.e. in the case when the notion of completeness 
does not apply. 

Completeness is not necessary for establishing the fact 
of correlations: if statistically significant correlations are 
found with an incomplete catalogue, this is a real signal. 
Incompleteness of the catalogue of BL Lacs cannot be a 
source of correlations with UHECR. Indeed, the objects 
which have to be added (removed) from a catalogue to 
make it complete are absent from (present in) the cata- 
logue for reasons not related to cosmic rays. Therefore, 
correlations with UHECR can only be weakened by such 
an incompleteness. 



B. UHECR autocorrelations 

The UHECR set used in our calculations is known 
to contain event clusters — in fact, it was chosen in 
such a way that auto-correlations are maximum. When 
studying cross-correlations of such a set with potential 
sources, one has to be careful to take clusters into account 
when calculating the probability of chance coincidence. 
Namely, one has to make sure that the signal observed 
in the data is not due to chance coincidence of clusters 
of cosmic ray events with candidate sources. For a given 
cluster, such a coincidence is roughly as probable as for a 
given single event, but contributes more into correlation 
function. This could produce artificial enhancement of 
correlations if not taken into account in the Monte-Carlo 
simulation. 

In our calculations |^ this problem is solved by in- 
troducing in each Monte-Carlo cosmic ray set the same 
number of doublets and triplets as there are in the real 
data. Then chance coincidences between clusters and 
candidate sources happen in simulated sets as often as 
in the real data, and thus are correctly accounted for. 
This is confirmed by calculations presented in Sect. Ill Dl 
where correlations with GRB catalogs do not appear de- 
spite auto-correlations in the UHECR dataset. 



C. Choice of angular scale 

The observed angular size of clusters suggests the an- 
gular scale at which correlations with sources are ex- 
pected. In our analysis Q we have fixed it to the previ- 
ously published HQ AGASA value of 2.5°. This value 
is treated as not adjustable. 

When the angular resolution of the experiment is 
known, there exists a preferred choice of the angular 
scale — the scale at which correlations are expected to 
be maximum. This angle can be determined by means of 
the Monte-Carlo simulation as follows. One has to gen- 
erate cosmic ray sets which are correlated with sources 
(taking into account the experimental angular resolution) 
and then measure angle at which the correlation signal is 
maximum. In the case of AGASA events this procedure 
gives the result very close to 2.5° (see dotted curve on 
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Fig. 3 of Ref. @), so in fact the choice of Ref. |^ was 
close to optimal. 

Alternatively, one may choose not to fix the angular 
scale and treat it as a free parameter. Then one should 
adjust it to maximize the correlation and calculate cor- 
responding penalty factor. A serious problem of this 
approach is that the result would depend on the lim- 
its within which the angular scale is varied. So finally 
one would have to input in one or the other way the in- 
formation about expected angular scale of correlations in 
order to obtain a definite answer. We do not follow this 
approach in our calculations. 



D. Hidden penalty 

We believe that we have accounted for all contribu- 
tions to the penalty factor ( "effective number of tries" ) 
directly related to our work. Still, the obtained signifi- 
cance may be somewhat overestimated. The point is that 
we are not the first who are looking for correlations be- 
tween UHECR and astrophysical objects. Some of these 
attempts have been published in the literature; the oth- 
ers may have never been reported. All these attempts 
should contribute, in principle, into the penalty factor. 
However, it does not seem possible to account correctly 
for all such contributions. 

The way around this problem is obvious. First, only 
very low values of P should be interpreted as a signal (in 
our calculations we considered P < 10""' to be sufficiently 
low to report our results). Second, and more important, 
the results have to be confirmed with a new independent 
data set. 



IV. FURTHER EVIDENCE 

The correlations found in Ref. for a particular set 
of UHECR and BL Lacs suggest, strictly speaking, only 
that UHECR and BL Lacs are connected to each other. 
They say little about the acceleration mechanism, parti- 
cle nature and other relevant physical parameters. These 
questions can be fully addressed only with much larger 
dataset than available now. Some attempts, however, 
can be made. The final purpose is to arrive at the end 
at a consistent picture which incorporates all known fea- 
tures of UHECR. Here is a sketch of these attempts, each 
bringing additional evidence of the connection between 
UHECR and BL Lacs. 

In Ref. (121 it was noted that cuts in BL Lac catalogue, 
chosen in Ref. |^ so as to maximize correlations with 
UHECR, select automatically 7-ray loud BL Lacs. When 
this observation is consistently elaborated, i.e. a sub- 
sample of the catalogue 13 is selected on the basis of 
a single criterium, namely, the cross-correlation with the 
EGRET sources, the resulting subset of 14 7-ray loud 
BL Lacs correlates with UHECR at the level of 10~^ 
of chance coincidence. This number cannot, of course. 



be interpreted as the significance of correlations between 
UHECR and BL Lacs because of a posteriori selection; 
rather the conclusion is that 7-ray loudness may be a 
distinctive feature of those BL Lacs which are UHECR 
accelerators. 

In Ref. an attempt was made to determine the 
charge composition of UHECR by reconstructing actual 
arrival directions of UHECR particles bent in the Galac- 
tic magnetic field. The idea was that such a reconstruc- 
tion should improve correlations of UHECR with BL Lacs 
if the latter are the sources. Substantial improvement 
was indeed observed for particles with the charge -1-1, 
which is an indication of the presence of protons. 

In Ref. |3| it was observed that if correlations of 
UHECR and BL Lacs were due to chance coincidence, 
the coinciding rays would be distributed over the sky ran- 
domly, reflecting only the local density of BL Lacs and 
exposure of a cosmic ray experiment. Thus, any signifi- 
cant deviation in the distribution of correlating rays over 
the sky from this expectation speaks in favor of real phys- 
ical connection between cosmic rays and BL Lacs. In fact, 
the UHECRs correlating with BL Lacs form two "spots" , 
with low probability to occur by chance |13|. This non- 
uniformity of the distribution of correlating rays may 
be due to several factors: (1) anisotropy of extragalac- 
tic magnetic fields at scales of order 500 Mpc; (2) poor 
knowledge of the Galactic magnetic field in some areas of 
the sky; (3) fluctuations in the space distribution of the 
nearest sources. 



V. CONCLUSIONS 

According to textbooks, any successful statistical anal- 
ysis consists in formulation of a "null hypothesis" and its 
subsequent falsification, at some confidence level, by com- 
paring to the experimental data. In the case at hand the 
"null hypothesis" which is being tested is that BL Lacs 
(and, in particular, any subset of them) and UHECR are 
uncorrelated. It takes only one counter-example to dis- 
prove a hypothesis, while "pro-examples" do not prove its 
validity. An illustration of this general rule has been dis- 
cussed in Sect. II: the absence of significant correlations 
of UHECR with the whole BL Lac catalogue does not 
prove that UHECR and BL Lacs are uncorrelated (i.e., 
that there is no subset of BL Lacs which are sources of 
UHECR and thus correlate with them). 

Having found such counter-example (i.e. the case when 
correlation is significant), the only thing which can be 
concluded, to a certain confidence level, is that UHECR 
and a particular subset of BL Lacs are correlated. The 
nature and physical implications of these correlations 
have to be studied separately by formulating and test- 
ing different "null hypothesis" . Refs.fTTlIT^ ll^ are first 
attempts in this direction. 
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APPENDIX A: STATISTICS OF CLUSTERING 
AND THE NUMBER OF UHECR SOURCES. 

The authors of Ref. [l| have misinterpreted our earlier 
paper J^] which concerns the statistics of clustering of 
UHECR and the number of their sources. In view of the 
growing confusion we would like to clarify this issue. 

In Ref. P] one reads: "(observed occurrence of clus- 
ters of UHECR) . . . was used to estimate the spatial 
density of sources to be 6 x 10~^ Mpc~^ 0. This 
would obviously place stringent constraints on candi- 
date astrophysical sources, e.g. 7-ray bursts (GRBs) 
have a spatial density of only ~ 10~^ Mpc~'^. However, 
a more careful analysis shows that the uncertain- 



ties in this estimate are very large. The true number 
is n = 2.77+2^53[2^7^) X 10-3 Mpc-3 at the 68% (95%) 
c.L; moreover relaxing the assumptions made, viz. that 
the sources all have the same luminosity and a spectrum 
cx E~'^, increases the allowed range even further, e.g. to 
n = 180™!?^;^) X 10-3 Mpc-3...". 

The actual situation is different from what is described 
in this extract. First, intrinsic inaccuracy of the estimate 
was fully acknowledged in Ref. Its physical reason 

is obvious: clustering is not sensitive to the number of 
dim sources. The estimates of Ref. cited above are a 
good illustration of the point: upper limits are huge and 
strongly model-dependent. 

Second, the main point of Ref. was to show that 
there exists a model-independent lower bound on the 
number of sources. This bound is presented in Ta- 
ble 1 of Ref. 0: at 90% and 99% confidence levels 
the number of sources is 71 > 2.3 x 10^"' Mpc-^ and 
n > 3.2 X 10^^ Mpc-3^ respectively. This is consistent 
with the calculations of Ref. [19j performed in particular 
models. 
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